Electronic polarity of nanoclusters: quantum and many-body effects 
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Interesting electrical polarity in nanoclusters usually requires the polarizability to exceed the value 
R 3 of the classical sphere of radius R. We clarify how this occurs naturally in single electron quantum 
systems, and relate it to the giant polarizability of NanFis, and to spontaneous dipole formation 
on niobium clusters. Many-body effects generally reduce the polarizability through screening. The 
usual RPA treatment retrieves the classical answer, but it significantly overestimates screening in 
few-electron systems. The system of two electrons on the surface of a sphere is solved numerically, 
to account for the Coulomb repulsion. At high densities, numerical results agree with RPA model 
with properly subtracted self-interaction effects. At low densities, the system performs quantum 
oscillations around the classical ground state. We calculate the lowest anharmonic correction to the 
polarizability, which also agrees well with numerical evaluation of the polarizability. 

PACS numbers: 73.22.-f, 77.80.-e 



I 

-a 
c 

o 
o 



> 

(N 

On 
O 

in 

o 

-i— > 
a 



i 

C 

O 

o 



Introduction. Electrical polaritji^ (both permanent 
moments /io and induced moments /Xi n d = <xF ex t) are 
important in determining the interaction of a nanoclus- 
ter with its environment^ and with external probes 
such as lights. Large polarizability may also lead 
to spontaneous symmetry-breaking, analogous to bulk 
ferroelectricity^Sii. Despite good numerical progress on 
the polarizability^ of small systems, its general under- 
standing is not developed. In classical electrostatics, a 
metallic sphere of radius R placed in an external electric 
field F develops surface charge whose field exactly can- 
cels the external field in the interior. The dipole fx = aF 
associated with this surface charge gives the polarizabil- 
ity of a classical metal sphere, «cms = R 3 - Experiments 
on various metallic clusters^ and on Cgcp^ yield values 
of a in rough agreements— However, larger a is neces- 
sary for spontaneous dipole formation. For example, in a 
crystal of nearly touching point-polarizablc spheres, the 
Clausius-Mossotti polarization catastrophe occurs when 
Airna/3 > 1, where n is the concentration. This trans- 
lates to a critical polarizability of a c = 3^/2R 3 /tt, larger 
than the available metallic value. 

How could larger-than-classical values of a emerge in 
metallic clusters? An analog case is the prediction^ and 
apparent confirmation^ that the system Na n +iCl n (and 
especially the highly symmetric Nai4Cli3 case), with 
one electron in a loosely-bound "surface state" outside 
a closed shell, has a second-order Jahn- Teller instability 
to a polar state. This suggests that the total polarizabil- 
ity (a i = electronic plus ^2 Z 2 s /uj 2 = vibrational, with 
Z e g being appropriate Born effective charges and lo the 
relevant vibrational frequencies^) has diverged. Simi- 
larly, de Heer's groups found low T permanent dipole 
moments on Nb clusters which might be interpreted as 
a(T) increasing to a divergence as T decreases. 

In this Letter, we analyze the polarizability of a nan- 
ocluster with a surface state occupied by an electron. As 



a model for this system, we consider an electron confined 
to the surface of a sphere. We show that due to quantum 
effects, electronic polarizability can considerably exceed 
the metallic value cxcms at low temperatures. This re- 
sult implies that the coupling between the surface elec- 
tron and elastic modes of nuclear displacement can de- 
crease the vibrational frequencies and drive the total po- 
larizability to diverge. We also analyze how the classical 
value acMS is recovered by screening when the number 
of metallic electrons increases. In particular, we consider 
the case of two surface electrons which is realized in the 
neutral cluster Mgi40i3. We predict that when a second 
electron is added, the polarizability diminishes to values 
closer to ckcms, and stability of the symmetric structure 
is restored. 

Polarizability of one electron on the sphere. Confin- 
ing electrons to the surface of a sphere permits theoret- 
ical simplification 15 . This model is not as much of a 
"spherical cow"— as might be thought. For example, the 
Ceo molecule confines 60 carbon tt electrons to the re- 
gion near the surface of a sphere. In neutral undistorted 
Nai4Cli3, one outer electron lies in a shell similar to a 
sphere. "Core-shell" nanoparticles with insulating cores 
and metallic shells have been studied— The electron 
eigenstates are ipL — Yl(SI), the spherical harmonics, 
with L = (£, m), and energies 
The polarizability a = (fi z ) JF Z is 
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where the dipole operator fi z is —eRcos 9 and wl — we — 
cxp(—Ei J /kBT)/Z is the probability that the system is 
in state L. At T = 0, the electron is in the ground 
state (wl = 5l.o)i an d there is only one non-zero matrix 
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FIG. 1: Polarizability versus inverse temperature in K for 
Nai4Fi3 clusters from Ref. ITsl . compared with theory for a 
single electron on a sphere. 



element, (1, 0| cos0|O, 0) = l/\/3. Then we get 



a(T = 0) = ^ 



(2) 



where as is the Bohr radius. This result is interesting, 
because it shows that for unscreened one-electron sys- 
tems, the value of a can exceed the value acMS = R 3 - 
The scaling with R 4 /as is not restricted to the spheri- 
cal shell, but is general for finite system with one elec- 
tron, e.g., a cubic box. The high temperature (ksT 3> 
h 2 /mR 2 ) answer can be found from classical statistical 
mechanics, and is given by Debye-Langevin law24: 



a (high T) 



e 2 R 2 
3k B T' 



(3) 



The crossover from the T = quantum answer to the 
classical answer occurs at crossover temperature fcgT t « 
H 2 /2mR 2 . Consider the outer electron of Nai4Fi3, re- 
sponsible for the giant polarizability measured by Rayane 
et ali&. The T-dependence can be interpreted as a large 
low T value from Eq. $2$, which evolves at higher T 
toward the classical value of Eq. @. At intermediate 
temperatures, the polarizability can be easily tabulated 
numerically. Using the well-known matrix elements of a 
dipole momenlii^ 



(I + 1, m\ cos8\l, m) = 
one can rewrite Eq. Q as 

t(T) = 6a(0)£ 
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The sum over m can now be done, giving 

a(T) = a(0) J2( w e - w w) = a (°) w o = a(Q)/Z. (6) 

e 

The computation reduces to calculation of the partition 
function Z = J2e( 2i + l)e~ £(£+1)/r , with the dimension- 
less temperature r = T/T*. Figure 1 shows values of 
a(T) from Eq. © plotted versus 1/r and compared with 
the dataiS for Nai4Fi3. Two fitting parameters were 
used: the value a(0) was taken to be 1950 A 3 , and the 
crossover temperature T* = h 2 /2mR 2 — 975K. Using 
Eq. (J2J, one finds reasonable values, R = 4.1 A and an 
effective mass m* — 5.3 free electron masses. (Of course, 
rotational and vibrational sources of T-dependence of a 
are probably not negligible, but the limited experimental 
range of T does not warrant a closer fit.) For compari- 
son, based on the Na-F bond length a = 2.3A, one can 
estimate the cluster size to be about i?=4A. Note that 
the low temperature polarizability is much larger than 
&CMS = 70 A 3 . 

One can consider an alternative explanation of the 1/T 
behaviour of the susceptibility. One can imagine that the 
cluster in fact has a large spontaneous dipole moment 
with fixed orientation with respect to cluster. Thermal 
rotations in the absense of an external electric field would 
randomize the direction of polarization, and an external 
field would align the dipole. Any mechanism of this sort, 
however, would produce saturation temperature much 
larger than T* because of a large mass associated with 
dipole rotations. This would be clearly incompatible with 
saturation temperature observed in Ref. Il9t 

In our calculation, we assumed that the response is lin- 
ear in the applied field. For this assumption to be valid, 
the polarization energy a£ 2 /2 must be smaller than the 
level spacing fi 2 /mR 2 . For larger fields, the dipole mo- 
ment would saturate. One can estimate the saturation 
field as f sa t = £ a t(as/i?) 3 , where £ at — e/a 2 B is atomic 
electric field, which is of the order of 10 10 V/m. One can 
see from this estimate that although the saturation field 
quickly decays with the cluster radius, it is still rather 
large for reasonable values of R. 

Screening effects. If there are many electrons on the 
shell, a mean field (RPA or time-dependent Hartree) ap- 
proximation may be reasonable. InA£*2°, polarizability 
was derived via RPA. Below we give a simple derivation 
and discuss the effects due to finite number of electrons. 
Consider the metallic sphere containing N electrons in 
the external field £ . The electric potential on the sur- 
face of the sphere is generated both by the external elec- 
tric field source and electrons. In polar coordinates, the 
external field contribution is — £r cos 6, while the dipole 
moment contribution is [i z cos9/r 2 for r > R. In the 
mean field theory, electrons respond to the screened po- 
tential on the surface r = R, which is proportional to 
the effective field £ c ff — £ — Hz/R 3 '■ The self-consistency 
condition is then 



V-z = a £ c 
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FIG. 2: Polarizability of two electrons on a sphere vs. its 
radius R, normalized to the classical value R 3 on a logarithmic 
scale. The mean field (Eq. © and semiclassical (Eq. t'Sll ) 
results are also shown for comparison. 



Solving Eq. for fi z , one arrives at 



<*MF 
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This calculation overestimates screening effects: the ac- 
tual field acting on an electron should not include the 
field created by the electron itself. Assuming that all 
electrons contribute equally to the screening field, one 
can exclude self-interaction by multiplying the screening 
field -ii 7 jR? by (1 - 1/iV). The result 



ttMFC — 



i + (i - 1/AOo-oAR 3 ' 



(9) 



unlike Eq. (JSJ, behaves correctly at N = 1. 

A realization of a dense metallic spherical shell is the 
7r-electron system of the Ceo molecule^. Pederson and 
Quong 21 computed the polarizability of Ceo both with- 
out and with self-consistent screening. Their answer for 
ao is 35 a? B per C atom, while after screening they find 
cxmf to be 9.3 a B per C atom. The electrons are on a 
spherical shell which is more than infinitesimally thin, 
but let us use the thin model anyway. To account for 
the reduction by screening, we need R — 9.0as. This 
compares reasonably with the radius 6.7-7.0 as of the 
carbon nuclei, adding an extra Angstrom of distance for 
the finite extent of the tt electron system. 

Two- electron model. The Coulomb interaction be- 
tween electrons leads to correlation between their posi- 
tions, which is omitted in mean field theory. To clarify 
the role of correlation, we consider two electrons on the 
sphere. Related studies of few-electron model systems 
but polarizability was not analyzed. We use 
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the Hamiltonian 
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n 2 



2mR 2 
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R ni - n 2 | 



(10) 



where ni i2 are positions of the electrons on the sphere 
{n\ 2 — 1), and V 2 2 are spherical Laplace operators. The 
first term in Eq. (|lf)|> represents the kinetic energy, and 
the second term describes Coulomb repulsion between 
electrons. The polarizability at T = is again given 
by Eq. Q}, with \L) being exact two-particle states, El 
their energies, and fi z — eR(n z \ + n z2 ) being total dipolc 
moment. There are multiple states with a given total 
angular momentum L' . The sum in Eq. Q is taken over 
all such states. 

Solving for eigenstates of the Hamiltonian l|10(l is dif- 
ficult even at zero temperature because of a vast con- 
figuration space. To reduce its dimension, we use sym- 
metry considerations. The ground state \l/o( n i) n 2) is 
symmetric with respect to electron coordinates and has 
total angular momentum £ — 0. Thus, the mass cen- 
ter of the two-electron system is completely delocalized, 
and ^o only depends on the distance between the elec- 
trons. We use the spherical angle 7 to parametrize this 
distance (cos 7 = ni • n 2 ), so the wave function takes the 
form \I/o( n i) n 2) = ^0(7)- To compute the action of the 
kinetic energy term V 2 2 in Eq. QlOfl on ^E'o( n i, n 2), it 
is convenient to use the following trick. Note that the 
kinetic energy in the £ = state also depends on the rel- 
ative angle 7 only, and therefore one can compute it at 
some particular values of ni, 2 and extend the result for 
all n 12 with the same value of 7. A convenient point cor- 
responds to the second electron at the north pole. Then, 
the angle 7 equals the polar angle of the first electron, 
and the kinetic energy is given by the well-known expres- 
sion for the Laplace operator on the sphere. Treating 
the second electron in the same manner, one obtains the 
Schrodinger equation for £ = states: 



Tt 1 d ( d \ 1 

m q- ( sin 7— ) + 0i=o(7) \ ¥0(7) = £0*0(7) , 

mn sin 7 07 \ 07 / J 

(11) 



where the effective potential 



C^=o(7) 



2i?sin^ 



(12) 



is entirely due to the Coulomb interaction. 

Now consider excited states. The only non-zero ma- 
trix elements in Eq. (JTJ are between I = and I — 1 
states. Moreover, since both ^0 and fJ-z are symmetric 
with respect to permutation of the electrons, one has to 
consider only those £ = 1 states in which electron mo- 
menta are added symmetrically. Up to a scalar multi- 
plier, one can construct only one symmetric vector out 
of ni j2 : N = (ni + n 2 )/(2cos^) (we normalized N so 
that N 2 = 1). In other words, the transformation prop- 
erties of £ = 1 states are identical to those of the mass 
center of the system. The wave function of the states 
with I = 1, £ z = therefore can be written in the form 

^1 (n.1,112) = N z ~&i (7). By a direct calculation, we 
find that the Schrodinger equation for ^1(7) has a form 
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similar to Eq. I|llf) with a different effective potential: 



C/ fc i(7) = U e =o{j) + 



AmR 2 



cos 2 3r 



(13) 



The last contribution in Eq. (|13fl is the kinetic energy of 
joint rotation of the two electrons with angular momen- 
tum L = 1. One can derive this "centripetal" term by a 
different method. At fixed distance between the electrons 
7, the rotational dynamics of the two-electron system can 
be described as that of a quantum asymmetric top with 
moments of inertia I\ = 2mR 2 sin 2 7/2, I2 = 2mR 2 , and 
I 3 = 2mi? 2 sin 2 7/2. It is known (see, e.g.^-) that i = 1 
is a special case in the quantum top problem: the wave- 
functions of these states are completely determined by 
the symmetry and therefore independent of -Zi.2,3- In our 
case it means that the top dynamics is independent of os- 
cillations of 7. This allows one to separate the variables 
for i = 1 state. The wave function of the two electron 
system can be written as \&i (111,112) = ^t(£, 0*7(7)) 
where is the wavefunction of the top described by 
Euler angles £, rj and £. (One can also check that this 
top wavefunction is proportional to N.) For symmetric 
1=1 states with zero projection of J on the top axis, 
the energy of the top is ^-(1/Ji + I//2), which coincides 
with the last term in Eq. Q13[l. 

One can now use Eqs. (|11|I - H13|I to find eigenfunctions 
and eigenvalues by solving a one-dimensional boundary- 
value problem. (The wave function has to be regular 
everywhere on the sphere. However, a general solution 
to Eq. may have singularities at 7 = and 7 = it, 
where the coefficients of Schrodinger equation are singu- 
lar. Thus, the requirement of regularity at these points 
serves as two independent boundary conditions.) While 
in the limit of R <C as one can thus reproduce the mean 
field theory result the opposite limit of strong in- 
teraction (R 3> as) is more interesting. In this regime, 
Coulomb repulsion dominates over kinetic energy. In the 
ground state, electrons perform zero-point oscillations at 
opposite poles of the sphere, thus minimizing potential 
energy. The electric field shifts the minimum and changes 
the energy of the system. The extra energy due to the 
field is t/g(7) = 2ef j_i?cos7/2, where £± is the field com- 
ponent perpendicular to the axis connecting electrons. 
Minimizing the potential energy L^ = o(7) + Us (7), one 
finds the ground state energy to be e 2 /2R — AR i £\. Av- 
eraging over all directions of the axis, one can replace £ 2 
by 2£ 2 /3. Thus, the polarizability of the two-electron 
system in the classical limit is 16i? 3 /3, still higher than 
the metallic value. 

However, it turns out that the amplitude of zero 
point motion around the equilibrium positions decays 
only as (as/R) 1 ^ 4 , and the simple classical picture de- 
scribed above holds only for impractically large values 
R/cib > 10 4 . To improve the classical approximation, 
one can expand Eqs. (|ll|l - l|13|) in the vicinity of the clas- 
sical equilibrium, 7 w 7r. For the £ = state, Schrodinger 



equation takes the form: 

h 2 



mR 2 



(h a + Sh )^ (x) =Eo^o{x) , 



(14) 



where x = n — 7 is the deviation from equilibrium. The 
dimensionlcss Hamiltonian 



hn = - 



dx 2 



ld_ 

x dx 



R 
2a^ 



(15) 



describes harmonic oscillations near the equilibrium, and 

5 R A 



x 

6ho = 3 & 



3 ■ 128 a B 



(16) 



is the first anharmonic correction to ho. For I = 1, one 
finds a similar equation with harmonic and anharmonic 
parts given by 



hi = h + 



1 



Shi = Sho + 



1 



(17) 



If one ignores anharmonic parts, these equations can be 
satisfied by wave functions of a harmonic oscillator: 



ip (x) 



where 



x 



exp(-x 2 /2xl) 
xip {x) 



(4a s /i?) 1 / 4 



(18) 
(19) 



(20) 




is the amplitude of zero point oscillations. The corre- 
sponding eigenvalues are 



2R\ V R j 



E a = — 1 + J — 
2R\ V R 

(21) 

The first term e 2 /2R is the classical Coulomb energy, and 
the second term is due to quantum oscillations. 

The zero point motion amplitude xq can be consid- 
ered as a small anharmonicity parameter. We shall now 
compute anharmonic corrections to the polarizability, ex- 
panding over the powers of Xq. Treating anharmonic 
terms 5ho and Shi by standard perturbation theory, one 
can compute the corrections to the energy levels: 



2 4 

e Xr\ 



(22) 



SE = (^oN'oH-^S 

R 16 

sei = (^\shM = e ^-$. 

R 16 



Similarly, one can also find corrections to the wave func- 
tions 
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24 
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9 



1 

36 



(23) 
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where £ = x/xq. 

One can now use 8E i and (h/>o.i to calculate the cor- 
rection to the polarizability given by Eq. 0J. To do that, 
one has to take into account anharmonic corrections to 
the dipole moment operator: 



fj, = 2eRN cos ^ w 2ei?Na;oC I 1 - 



24 



(24) 



(here N is the unit vector directed to the mass center), 
and to the volume element in the integral: 



dQ oc sin jdj 



e4 

6 



adz 



(25) 



(One can check that this volume element is consistent 
with kinetic part of the Hamiltonian.) In the lowest or- 
der in anharmonic corrections, the only important con- 
tribution comes from the lowest 1=1 state. The con- 
tributions of higher I = 1 states are absent in harmonic 
approximation, and their contributions to the matrix el- 
ement are proportional to anharmonicity parameter x\. 
Since the polarization contains square of the matrix ele- 
ment, excited J = 1 states contribute only in the fourth 
order in xq, which we do not consider here. Substituting 
Eq. H24|) into Eq. and replacing the matrix element of 
N z between J = and J = 1 states by one finds: 



2e 2 R 2 



M 



3 Ei + 8 Ex -E - 8E NqNx 



(26) 



where 



M = 2 / 1 



24 



(^0 + ^0)^1+^1) (27) 



and iVo,i are normalization integrals, with corrected vol- 
ume element: 



v, = -> / (l-^)t<-,-.w,r< ( /< . (2M 



Substituting wave function corrections (|23|l . one finds, in 
the lowest order: 



(29) 



M = 1 _ H 
3 



The energy difference, according to Eqs. (|21|l and (|22|l . 

is 



£?! - So + «5i - SE = ^ C X l + M 



(30) 



Substituting Eqs. and (^UJ into Eq. I2BJ) , one finally 
obtains: 



o.i(R > as) 



16i? 3 
3 

16i? 3 



49 2 
36 Xn 



(31) 



49 /asX 1 / 2 
18 V~R ) 



The results of numerical solution of the two-electron 
problem are shown on Fig. [21 The self-consistent re- 
sult © breaks down at R ~ 3as, while the classical 
regime <|31[1 takes over at R > lOOa^. For reasonable 
size nanoclusters (R ~ lOa^), the effect of electron cor- 
relations is significant, but the polarizability is still con- 
siderably higher than the classical value a = i? 3 . 

The limit R < clb shown in the left part of the plot may 
seem unphysical. However, note that the Bohr radius 
entering our calculation is in fact an effective parameter, 
depending on the effective mass of the surface state. If, 
for some reasons, the effective state is a light electron, 
the effective Bohr radius can also be large, and the limit 
R < as may become physically applicable. 

It is also instructive to check whether the Clausius- 
Mossotti catastrophe (Air/ina > 1) can be achieved. 
Packing spheres in the most dense way, one requires 
a/R 3 > 3\f2/n, which is achieved for R > 2.5as. Thus, 
despite screening, an array of two-electron nanoclusters 
can develop a spontaneous polarization. 

To summarize, the polarizability of a nanocluster with 
a surface electron state may be significantly enhanced 
due to quantum effects. This enhancement, however, dis- 
appears when the number of "metallic" electrons is in- 
creased. Although we considered a very simple model of 
a nanocluster, our results should hold qualitatively for 
more complex clusters, including niobium ferroelectric 
clusters studied in Ref.|(| At the first sight, these clusters 
do not look like a single-electron systems. However, if the 
spectrum of the high T state has its outermost electron 
in a separated level with orbital degeneracy, as in Na3 
(see Ref fl4|l . or with quasi- Jahn- Teller degeneracy as in 
Nai4Cli3, similar physics probably applies. Also, despite 
the fact that we only considered surface electrons, one 
can expect a similar physics when electrons are confined 
in the bulk of the cluster. For example, the estimate for 
the level spacing, dipole moment and therefore for polar- 
izability, should be still valid for bulk electrons. Thus, 
our results can give correct qualitative description for 
more complex clusters with a small number of "mobile" 
electrons. 
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